Quantum Phase Diffusion in a Small Underdamped Josephson Junction 
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Quantum phase diffusion in a small underdamped Nb/A10_v/Nb junction (~ 0.4 yum^) is demonstrated in a 
wide temperature range of 25-140 mK where macroscopic quantum tunneling (MQT) is the dominant escape 
mechanism. We propose a two-step transition model to describe the switching process in which the escape rate 
out of the potential well and the transition rate from phase diffusion to the running state are considered. The 
transition rate extracted from the experimental switching current distribution follows the predicted Arrhenius 
law in the thermal regime but is greatly enhanced when MQT becomes dominant. 

PACS numbers: 74.50.+r, 05.40.-a, 85.25.Cp 
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Classical and quantum diffusion of Brownian particles in 
titled periodic potential plays a fundamental role in the dy- 
namical behavior of many systems in science and engineering 
1 1-16]. Examples include current biased Josephson junctions 
I jj-lsd, colloidal particles in arrays of laser traps iflollTlll. cold 
atoms in optical lattice or Bose-Einstein condensates 112141411 . 



and various biology-inspired systems known as Brownian mo- 
tors (molecular motors or life engines), which receive consid- 
erable attention in physics f\3\ and chemistry fl^ . Because 
of the design flexibility, manufacturability, and controllability 
Josephson junctions provide an excellent test bed for making 
quantitative comparison of experimental data with theoretical 
predictions and unraveling possible new physics in the tilted 
periodic potential systems. 

The dynamics of a cuiTent biased Josephson junction can 
be visualized as a fictitious phase particle of mass C moving 
in a tilted periodic potential U{ip) - -Ej(i(p + cos (p). Here, C 
is junction capacitance, / = I/Ic is the junction's bias current 
normalized to its critical current, the phase particle's position 
(p is the gauge invariant phase difference across the junction, 
and Ej = Tildle is the Josephson coupling energy with e and fi 
being the electron charge and Planck's constant, respectively. 
Previous experiments using Josephson junctions have identi- 
fied three distinctive dynamical states, as shown schematically 
in Fig. 1. In the first state, the phase particle is trapped in one 
of the metastable potential wells and undergoes small oscil- 
lation around the bottom of the well with plasma frequency 
(Op. Because of thermal and/or quantum fluctuations the par- 
ticle has a finite rate Fi escaping from the trapped state. The 
escape rate becomes significant when the barrier height At/ 
is not much greater than ksT or fitOp, where kg is the Boltz- 
mann constant and T denotes the temperature, respectively. 
After the particle escapes from the initial well, depending on 
the energy gain 6U - ('I'o being the flux quantum) and 
the loss Eo due to damping (cf. Fig. 1), it could enter either 
the second dynamical state called phase diffusion (PD) or the 
final running state. In the former case as the bias cuiTent / 
is increased further the particle will eventually make a transi- 
tion, characterized by a rate constant F2, to the running state. 



While escape from the trapped state to PD is difficult to detect 
transition to the running state is signaled by a sudden jump in 
the dc voltage of the junction (called switching) and thus can 
be readily captured in real time by increasing / continuously 
from zero until a switching occurs iITtIi . 

The fundamental importance of understanding PD has stim- 
ulated many studies in recent years. However, experimental 
studies were focused mostly on the classical regime where 
thermal activation (TA) is the dominant escape mechanism 
and thermal fluctuation governs the PD process On the 

other hand, in the quantum regime where macroscopic quan- 
tum tunneling (MQT) dominates, one expects that quantum 
fluctuation induced tunneling will play an important role in 
the PD process and subsequent transition to the running state 
thus the term quantum PD (QPD) has been coined in the lit- 
erature Iil5t, il8t-t20il . However, although theoretical progress 
of QPD in overdamped systems has been remarkable over re- 
cent years the situation is so far much less clear for 
underdamped systems [15, 20]. 

In this work, we demonstrate QPD in a small underdamped 
Josephson junction over a wide temperature range of 25 to 140 
mK. To contrast QPD with classical PD, we use two Nb-AlO j- 
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FIG. 1: (Color online) Phase particle in the trapped, diffusion, and 
running states (denoted by n = 1,2, 3, respectively) with occupation 
probability p„ in a tilted washboard potential. 
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Nb trilayer junctions of different sizes (see Table I) having Tq 
«: Tcr and Tq » T^r, respectively. Here, To is the temperature 
above which PD occurs and Tc- is the classical-to-quantum 
crossover temperature below which MQT dominates. One of 
the hallmarks of PD in underdamped junctions is the narrow- 
ing of the width cr of switching current distribution P(I) as 
temperature increases ||5|-l8|]. This is observed clearly in the 
measured cr{T) of the larger junction L above Tq ^ 450 mK 
» T^i; which indicates that PD in this case is classical in na- 
ture. In sharp contrast, for the smaller junction S the width cr 
continues to increase as temperature decreases to the lowest 
value of 25 mK. When plotted in semilogarithmic scale cr vs 
T shows a clear increase of slope around T^^. = 140 mK, point- 
ing to a change from classical PD to QPD. We will extract the 
transition rate r2 directly from the experimental results and 
show that QPD is fundamentally different from classical PD. 

Two Nb/A10j:/Nb junctions used in this study were fab- 
ricated on the same chip with nominal areas of 0.52 and 
1.61 yum^ for junctions S and L, respectively. Compared 
with previous works reported in Refs. (3] and |@], where dc 
SQUIDs were used to tune Ic, our approach kept Ic/C con- 
stant. This unique approach is essential to extend PD to 
the quantum regime. Since r„- = fta)p[(l + 1/40^)'^^ - 
I /2Q]/27TkB ~ hdJo/lnkB scales with the plasma frequency 

(jjp - (jjo (l - f^^''^ , where ojq - {Inld^oC)^'^ and Q - 
(jjpRC {R being junction's damping resistance), Tcr is approx- 
imately independent of the junction sizes as long as they are 
fabricated from the same trilayer On the other hand, To can be 
reduced by making smaller junctions therefore we are able to 
tune To and Tcr independentl y to meet the condition Tq <sc Tcr 
required for observing QPD 112 111 . 

Figure 2 shows the measured P{I) from 25 to 800 mK for 
junction S with its I-V curve at 30 mK displayed in the in- 
set. In our experiment, P{I) was measured by the time-of- 
flight technique ^Jt^ with dildt - 1 10/sec for sample S and 
163/sec for sample L. Each measured P{I) consisted of 50000 
switching events. In Fig. 3, we plot cr and the mean Is of P(I) 
versus temperature (symbols) for junction S in (a) together 
with those of junction L in (b). For junction L the measured 
cr{T) shows the famihar classical PD started at temperature T^ 
^ 450 mK well above T^^^ = 125 mK. The solid lines in (b) are 



TABLE I: Parameters of two Nb/AlO j/Nb junctions S and L used in 
this work. Rn is normal-state resistance obtained from /-V curves. I,., 
C, and R for L are determined from fits to experiment using TA and 
MQT theories below 450 mK and Monte Carlo simulations above it. 
Those for S are obtained considering its Rn ratio to L (Note a slightly 
larger R chosen to have a better fit). See the text for details. 

Junction Area''0"m') Rn{U1) /,(nA) C(fF) R{n) r„(mK) ro(mK) 



S 
L 



0.39 
1.54 



15.1 
3.84 



122 
480 



19.6 1800 
77 315 



140 
125 



<25 
-450 



"Estimated for L from fitted C and a specific capacitance of 50 fF//jm'. The 
value for S is obtained via its Rfj ratio to L. Nominal areas for junctions S and 
L were 0.52 and 1.61 fim^, respectively. 
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FIG. 2: (Color online) Experimentally measured P(I) of junction S 
at some temperatures indicated. The inset shows the I-V trace of the 
junction at 30 mK. 



calculated according to the TA 112311 and MQT 112411 rate formu- 
las using the parameters listed in Table I. The dashed lines are 
from Monte Carlo simulations considering thermal fluctuation 
and PD isK^j- In contrast to junction L the observed cr for 
junction S in Fig. 3(a) shows a monotonic decrease with in- 
creasing temperature, indicating that PD occurred in the entire 
temperature range of the experiment. Furthermore when plot- 
ting the data in semilogarithmic scale as shown in the inset of 
Fig. 3 we notice a distinctive slope decrease around T;?^ = 140 
mK from MQT to TA regimes. Such a decrease can be eas- 
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FIG. 3: (Color online) (a) Width cr and mean /j of experimental P{I) 
of junction S (symbols), (b) Corresponding data of junction L. Solid 
lines in (b) are calculated from TA and MQT theories while dashed 
lines from Monte Carlo simulations considering thermal PD |5, 8]. 
The inset shows cr of junction S plotted in semilogarithmic scale. 
Two solid lines are guides to the eye displaying a slope turning near 
Ti = 140 mK. 
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ily understood since TA causes cr to increase with increasing 
T which partially cancels the effect of negative (l / o-)dcr / dT 
due to PD. 

To gain further insight and have a quantitative grasp on the 
effects of escape (from the trapped state to PD) and transition 
(from PD to the running state) on switching current distribu- 
tion, regardless of whether TA or MQT is the dominant mech- 
anism, we set up the following master equation according to 
the two-step transition model shown in Fig. 1: 



dpi/dt - 
dp2/dt — 
dp^/dt — 



= -Ti pi 

= Ti pi - Fa p2 

= P2, 



(1) 



where p„ (n - 1, 2, 3) is the probability of finding the phase 
particle in state n. Since P(/) = dp^/dl, it follows straightfor- 
wardly that 



r2(/) 



(dI/dt)P(I) 



(2) 



Equation (|2]) shows that F2(/) can be extracted from measured 
P(I) provided Fi (/) is known, which is true in our experiment. 
Notice that in the limit of r2 — ^ 00, Eq. (|2]i leads directly to 
Fi(/) = {dI/dt)P{I)/[l - £ P(l')dl'] which is identical to the 
result of Fulton and Dunkleberger IftIi in which PD is ab- 
sent. In the opposite limit of F2 Fi, the same expression 
is obtained with Fi replaced by F2: F2(/) = (dl /dt)P(I)/[l - 

rl I I 

Jq P{I )dl ]. These results mean that the much slower process 
plays the major role in determining P{I), as expected. In the 
more general situation of F2 ~ Fi, Eq. (2) enables one to sep- 
arate the effect of F2 on switching current distributions from 
that of Fi. The inverse procedure of computing P{I) from Fi 
and F2 is given by: 



P{I) 



(dl/dty 



Jo 



(3) 



Equations (|2]) and (O thus allow us to quantitatively investi- 
gate the dependence of (Q)PD on bias current and the inter- 
play between the particle's escape and (Q)PD. In Fig. 4(a), we 
plot Fi (solid lines) calculated using the parameters of junc- 
tion S and F2 (symbols) extracted from the measured P(/) us- 
ing Eq. It can be seen that at T = 800 mK, Fi is several 
orders of magnitude greater than F2. The measured P(/) is 
therefore entirely determined by F2. As the temperature de- 
creases, Fi is seen to progressively approach F2. 

Having clearly established that PD occurs in both classical 
and quantum regimes in junction S, we now use the data in 
Fig. 4(a) to further demonstrate the key difference between 
classical PD and QPD. In Fig. 4(b), we plot F2 versus l/T 
at three bias currents (thus fixed potentials) of 48, 52, and 56 
nA, which shows distinct features below and above 7;?^. While 
the data above T^,. follow the straight lines, indicating that F2 
in the classical regime obeys the Arrhenius law F2 displays a 
much weaker IjT dependence at below T <s: T^^. We note that 
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FIG. 4: (Color online) (a) Transition rate Y2 (symbols) and escape 
rate Fi (solid lines) of junction S at some typical temperatures, (b) 
r2 ~ 1/r at three fixed currents as indicated by the vertical arrows in 
(a). Dashed and dotted lines are fits displaying the Arrhenius law. (c) 
/ ~ r for fixed r2 = 2000 sec"' as indicated by a horizontal arrow in 
(a). Solid lines in (b) and (c) are guides to the eye. 



similar behavior in the classical regime was discussed previ- 
ously by Vion et al. [3] for overdamped system where the 
diffusive particle is considered to overcome an effective dissi- 
pation barrier In that case, the transition rate from PD to the 
running state, which retains the familiar Kramers form, was 
derived. Fitting the data above T^,. using F2 = a exp(-b/T), 
we obtain fl = 5.2x10^ sec"', = 2.3 Kfor / = 48 nA (dashed 
hne) and a = 3.3 x 10** sec"', ^ = 1.7 K for / = 52 nA (dot- 
ted line). The effective barrier b appears smaller as compared 
to the calculated barrier height AU of 2.68 and 2.46 K due to 
the motion of the diffusive particles, which is physically quite 
reasonable. These results indicate that in the thermal regime a 
dissipation-barrier description is also appUcable to PD in un- 
derdamped junctions. 

Machura et al. recently investigated the diffusion problem 
of overdamped particles using the Smoluchowski equation in- 
corporating quantum fluctuations |19]. They found that the 
particle's average velocity (v) increases with increasing tem- 
perature and quantum effects always assist the particle to over- 
come barriers leading to a larger (v) than that in absence of 
quantum fluctuations. Because in our underdamped junction 
the dc voltage, which is proportional to (v), produced by PD 
is too low to be detected directly |25], it can nevertheless be 
expected that a larger (v) would result in a larger F2 since the 
increased kinetic energy makes transitions to the running state 
easier. For this reason, the data in Fig. 4(b) are consistent with 
the theoretical prediction since extrapolating F2 from the clas- 
sical to the quantum regime would lead to rates that are much 
lower than the experimental data. Therefore, the much weaker 



4 



1 / r dependence of r2 below 7"^^, in a stark contrast to the Ar- 
rhenius behavior above T^^, manifests the quantum nature of 
the diffusion process atT < T^^. 

In Fig. 4(c) we plot / versus T for a constant r2 - 2000 
sec"\ which again shows a distinctive change of slope around 
similar to that of cr. The approximate linear / de- 
pendence above T^^ can be qualitatively explained. In the 
absence of thermal fluctuations transition from PD to run- 
ning state is expected to occur deterministicaUy at Iq where 
6Uq - (h/2e)Io - Ed- For 7" > the phase particle will 
exit the PD state prematurely because the particle on aver- 
age acquires an additional thermal energy of ~ ksT. Thus 
the condition for transition out of PD needs to be revised to 
6U + ksT = Ed. Assuming the junction's damping, and thus 
Ed, saturates at low T we obtain ih/2e)I = Ed - ksT. The 
predicted slope |s| - lekg/h » 7 nA/K is comparable to 
the experimental value of 15 nA/K in the thermal regime in 
Fig. 4(c), which is quite reasonable considering the simplicity 
of the model. Below T^^,, however, the measured |^| increased 
to about 68 nA/K, about an order of magnitude greater than 
leks/h which remains unexplained. 

In conclusion, QPD was demonstrated and systematically 
studied in a small underdamped Nb Josephson junction. Us- 
ing junctions of different sizes fabricated on the same chip 
we were able to calibrate the relevant parameters of the small 
junction and at the same time extended QPD over a wide tem- 
perature range. We showed that cr decreases monotonically 
with increasing temperature and there is a distinctive change 
of slope at Tcr below and above which QPD and classical PD 
occur. We developed a two-step transition model with which 
the effects of escape rate Fi (from the trapped state) and the 
transition rate (from PD to the running state) on switching 
current distributions can be separated and F2 be determined 
from the measured P(I) directly. It was found that F2 vs T at 
fixed bias current, and thus fixed potential landscape, follows 
the Arrhenius law in the case of classical PD. The most im- 
portant finding was that for QPD, r2 is exponentially higher 
than that expected for the classical PD and has a much weaker 
l/T dependence. The similarities between the temperature de- 
pendence of Fi and F2 in underdamped Josephson junctions 
going from classical regime to quantum regime were strik- 
ing. We hope our experimental progress and advancement in 
data analysis will stimulate further theoretical and experimen- 



tal studies of and lead to a better understanding of the quantum 
diffusion phenomena in underdamped tilted periodic potential 
systems. 
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